Theory and simulation of two-dimensional nematic and tetratic phases 
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Recent experiments and simulations have shown that two-dimensional systems can form tetratic 
phases with four-fold rotational symmetry, even if they are composed of particles with only two-fold 
symmetry. To understand this effect, we propose a model for the statistical mechanics of particles 
with almost four-fold symmetry, which is weakly broken down to two-fold. We introduce a coefficient 
K to characterize the symmetry breaking, and find that the tetratic phase can still exist even up to a 
substantial value of k. Through a Landau expansion of the free energy, we calculate the mean-field 
phase diagram, which is similar to the result of a previous hard-particle excluded-volume model. 
To verify our mean-field calculation, we develop a Monte Carlo simulation of spins on a triangular 
lattice. The results of the simulation agree very well with the Landau theory. 

PACS numbers: 64.70.mf, 61.30.Dk, 05.10.Ln 



I. INTRODUCTION 

In statistical mechanics, one key issue is how the mi- 
croscopic symmetry of particle shapes and interactions 
is related to the macroscopic symmetry of the phases. 
This issue is especially important for liquid-crystal sci- 
ence, where researchers control the orientational order 
of phases by synthesizing molecules with rod-like, disk- 
like, bent-core, or other shapes. In many cases, the low- 
temperature phase has the same symmetry as the parti- 
cles of which it is composed, while the high-temperature 
phase has a higher symmetry. For example, in two di- 
mensions (2D), particles with a rectangular or rod- like 
shape, which has two-fold rotational symmetry, form a 
low-temperature nematic phase, which also has two-fold 
symmetry. Likewise, if the particles are perfect squares, 
which have four-fold rotational symmetry, they can form 
a four-fold symmetric tetratic phase. 

An interesting question is what happens if the symme- 
try of the particles is slightly broken. Will the symmetry 
of the phase also be broken, or can the particles still form 
a higher-symmetry phase? For example, we can consider 
particles with approximate four-fold rotational symme- 
try that is slightly broken down to two-fold, as in Fig. [TJ 
Can these particles still form a tetratic phase, or will they 
only form a less symmetric nematic phase? 

Recently, several experimental and theoretical studies 
have addressed this problem. Narayan et al. [l| per- 
formed experiments on a vibrated-rod monolayer, and 
found that two-fold symmetric rods can form a four- 
fold symmetric tetratic phase over some range of pack- 
ing fraction and aspect ratio. Zhao et al. p] stud- 
ied experimentally the phase behavior of colloidal rect- 
angles, and found what they called an almost tetratic 
phase. Donev et al. Q simulated the phase behavior of 
a hard-rectangle system with an aspect ratio of 2, and 
showed they form a tetratic phase. Another simulation 



by Triplett et al. 3 showed similar results. In further 
theoretical work, Marti'nez-Raton et al. Q developed 
a density-functional theory to study the effect of particle 
geometry on phase transitions. They found a range of 
the phase diagram in which the tetratic phase can exist, 
as long as the shape is close enough to four-fold symmet- 
ric. In all of these studies, the particles interact through 
hard, Onsager-like excluded-volume interactions. 

The purpose of the current paper is to investigate 
whether the same phase behavior occurs for particles 
with longer-range, soft interactions. We consider a gen- 
eral four-fold symmetric interaction, which is slightly bro- 
ken down to two-fold symmetry. We first calculate the 
phase dia gram using a Maier-Saupe-like mean-field the- 
ory [qI. Il0||. To verify the theory, we then perform 
Monte Carlo simulations for the same interaction. 

This work leads to two main results. First, the tetratic 
phase still exists up to a surprisingly high value of the 
microscopic symmetry breaking (as characterized by the 
interaction parameter k, which is defined below). Sec- 
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FIG. 1: Schematic illustration of an interacting particle sys- 
tem in the tetratic phase. The shape of the particles indi- 
cates that the rotational symmetry of the interaction is bro- 
ken down from four-fold to two-fold. 
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ond, the phase diagram is quite similar to that found by 
Martinez- Raton et al. for particles with excluded- volume 
repulsion. This similarity indicates that the phase behav- 
ior is generic for particles with almost-four-fold symme- 
try, independent of the specific interparticlc interaction. 

The plan of this paper is as follows. In Section [TTl we 
present our model and calculate the mean-field free en- 
ergy. We then examine the phase behavior and calculate 
the phase diagram in Section Hill In Section |TV] we de- 
scribe the Monte Carlo simulation methods and results. 
Finally, in Section |V] we discuss and summarize the con- 
clusions of this study. 

As an aside, we should mention one point of terminol- 
ogy. The tetratic phase has occasionally been called a 
"biaxial" phase, by analogy with 3D biaxial nematic liq- 
uid crystals (5|. However, this analogy is somewhat mis- 
leading. In 3D liquid crystals, the word "biaxial" refers 
to a phase with orientational order in the long molecular 
axis and in the transverse axes, i.e. a phase with lower 
symmetry than a conventional uniaxial nematic. By con- 
trast, the tetratic phase has higher symmetry than a con- 
ventional nematic, four-fold rather than two-fold. For 
that reason, we will not use the term "biaxial" in this 
work. 



II. MODEL 

Maier-Saupe theory is a widely used form of mean-field 
theory, which describes the isotropic-nematic transition 
in 3D liquid crystals. In this section, we extend Maier- 
Saupe theory to describe 2D liquid crystals with almost- 
four-fold symmetry, as shown in Fig.[T] For this purpose, 
we use the modified Maier-Saupe interaction 

Ul2 (ri2, 012) = -C/o(ri2) [« cos (2012) + cos (4012)] , (1) 

where 0i2 = 0i — 02 is the relative orientation angle 
between particles 1 and 2, and ri2 is the distance be- 
tween these particles. In this interaction, the dominant 
orientation-dependent term is cos(40i2), which has per- 
fect four- fold symmetry. The term cos(20i2) represents 
a correction to the interaction, which has only two-fold 
symmetry. If the coefficient k is small, then the sym- 
metry is slightly broken from four-fold down to two-fold. 
(By contrast, if k is large, then the interaction clearly 
has two-fold symmetry and the four-fold term is unim- 
portant, as in classic Maier-Saupe theory.) The over- 
all coefficient C/o(?'i2) is an arbitrary distance-dependent 
term. 

In mean-field theory, we average the interaction energy 
to obtain an effective single-particle potential due to all 
the other particles, 



which is normalized as 



Po 



d0p(0). 



(3) 



f/off (0i) = / d2ri2d02P (02) Ui2 (ri2, 012) 



(2) 



Here, p (02) is the orientational distribution function, 



where po is the number density of particles. To calculate 
Ucs, we set the x-axis along an ordered direction (the 
director in nematic case, or one of the two orthogonal 
ordered directions in the tetratic case). In that case, the 
averages of sin(20) and sin(40) vanish by symmetry, and 
hence Eq. ^ becomes 

f/cff (0) = -Upo [kC2 cos (20) + Ci cos (40)] , (4) 

where U is the integral over the position-dependent part 
of the potential, and 

C2 = (cos (20)), (5a) 
C4 = (cos (40)). (5b) 

The resulting orientational distribution function is 
^ Pa exp[7(KC2 cos(20) + C4 cos(40))] 
^ ^ lode exp[7(KC2 cos(20) -I- C4 cos(40))] ' ^ ' 

where we have defined the dimensionless ratio 7 — 
poU/{kBT). 

Note that C2 can be regarded as a nematic order pa- 
rameter, and C4 as a tetratic order parameter. In the 
isotropic phase, the system has C2 = C4 = 0. By com- 
parison, in the tetratic phase, the system has C2 — but 
C4 7^ 0. In the nematic phase, with the most order, the 
system has C2 7^ and C4 7^ 0. 

To determine which of these phases is most stable, we 
must calculate the free energy F = (U) + kBT (log p) as a 
function of the order parameters C'2 and C4 . The average 
interaction energy per particle is 



{U) = -^Upo {nCl + Cl) 



(7) 



The entropic part of the free energy per particle is 
kBT (log p) = Upo {i^Ci + Cl) (8) 

-ksT log d0exp[7(KC2Cos20-|-C4COs40)]^ ; 

here we have have subtracted off the constant entropy of 
the isotropic phase. We combine these terms and nor- 
malize by fc^T to obtain the dimensionless free energy 

^ = i7K|H-CS (9) 

- log ^- ^ d0 exp[7(KC2 cos 20 + C4 cos 40)] 

Minimizing this free energy with respect to C2 and C4 
gives the equations 

C2 = — [ d0cos(20)p(0) , (10a) 
po Jo 

C4 = — [ d0cos(40)p(0), (10b) 
Po Jo 

which are exactly consistent with Eqs. 
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III. MEAN-FIELD RESULTS 

The model is now completely defined by two dimen- 
sionless parameters: 7 — poU / [kBT] is the ratio of in- 
teraction energy to temperature, and k represents the 
breaking of four-fold symmetry in the interparticle inter- 
action. We would like to determine the phase diagram 
in terms of these two parameters. As a first step, we 
minimize the free energy of Eq. ^ numerically using 
Mathematica. We then do analytic calculations to ob- 
tain exact values for second-order transitions and special 
points in the phase diagram. 

Figure [2] shows the numerical mean-field results for 
the order parameters C2 and C4 as functions of 7, for 
several values of k. These plots represent experiments 
in which the temperature is varied, for particles with 
a fixed interaction. When the four-fold symmetry is 
only slightly broken by the small value k = 0.4, there 
are two second-order transitions, first from the high- 
temperature isotropic phase to the intermediate tetratic 
phase, and then from the tetratic phase to the low- 
temperature nematic phase. For a larger value k = 0.75, 



the isotropic-tetratic transition is still second-order, but 
now the tetratic-nematic transition is first-order, with a 
discontinuous change in C2. For k ~ 1.5, the two tran- 
sitions merge into a single first-order transition directly 
from isotropic to nematic, with discontinuities in both C2 
and C4, and the tetratic phase does not occur. Finally, 
for the largest value k = 2.25, the isotropic-nematic tran- 
sition becomes second-order; this behavior corresponds 
to the prediction of 2D Maier-Saupe theory with a sim- 
ple cos 2^12 interaction. 

The numerical mean-field results are summarized in 
the phase diagram of Fig. [31 The system has an isotropic 
phase at low 7 (high temperature) and a nematic phase 
at high 7 (low temperature). It also has an intermedi- 
ate tetratic phase, as long as the symmetry-breaking k is 
sufficiently small. The temperature range of the tetratic 
phase is very large for small k, then it decreases as k in- 
creases, and finally vanishes at the triple point B. In this 
mean-field approximation, the isotropic-tetratic transi- 
tion is always second-order and independent of n. The 
tetratic-nematic transition is second-order for small k, 
then becomes first-order at the tricritical point A. The di- 
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FIG. 3: Phase diagram of the model in terms of 7 (inverse 
temperature) and ft (two- fold distortion in the interaction). 
The grey solid lines represent second-order transitions, and 
the dark solid lines are first-order transitions. The dashed 
lines indicate the extrapolated second-order transitions, which 
give the cooling limits of the metastable phases. Point B 
(0.79,2) is the triple point, and A (0.61,2.2) and D (2,1) are 
the two tricritical points. Point C (1,2) is the intersection of 
the extrapolated second-order transitions. 



rect isotropic-nematic transition is second-order for large 
K, then becomes first-order at the tricritical point D. 
Point C is the intersection of the extrapolated second- 
order transitions, and represents the limit of metastabil- 
ity of the tetratic phase. 

To calculate second-order transitions and special 
points in the phase diagram, we minimize the free energy 
of Eq. © analytically. For this calculation, we expand 
the free energy as a power series in the order parameters 
C2 and C4, which gives 



F 



7k(2-7k) 2 , 7(2-7) 



8 

44 4 



cl 



— ^L.204 



(11) 



Note that this expression is exactly what would be ex- 
pected in a Landau expansion based on symmetry; it is 
always an even function in C2, but it is an even function 
of C4 only when C2 = 0. 

To find the isotropic-tetratic transition, we set C2 = 
in the expansion, because this order parameter vanishes 
in both of those phases. The second-order isotropic- 
tetratic transition then occurs when the coefficient of C| 
passes through 0. Hence, the transition is at 



7 



(12) 



independent of n. 

For the second-order isotropic-nematic transition, we 
see that the isotropic phase becomes unstable when 
a^F/aCl = d'^F/dCl = d'^F/dC2dCi = 0, all evaluated 
at C2 = C4 = 0. These equations have two solutions. 



FIG. 4: A snapshot of the spins on a triangular lattice in the 
tetratic phase. The shape of the rectangles is just a schematic 
illustration of the symmetry of their interaction. 



one of which corresponds to the isotropic-tetratic tran- 
sition found above. The other solution, representing the 
isotropic-nematic transition, is 



(13) 



On the nematic side of this transition, we find C4 = 
K^7^C|/[4(2 — 7)]; i.e. the order parameters C2 and C4 
increase with different critical exponents. We substitute 
that relation into the expansion (|lip to obtain an effective 
free energy in terms of C2 alone. 



Faff 



7k(2 — 7k) 
8 



7"^ni-7) 
32(2-7) 



(14) 



The tricritical point D occurs when the coefficients of both 
C| and C| vanish in this expansion, which is at 7 = 1 
and K = 2. 

For the second-order tetratic-nematic transition, we 
cannot use the expansion of Eq. p4|) because C4 is not 
necessarily small; instead we return to the free energy of 
Eq. Anywhere in the tetratic phase we must have 
dF/dCi — 0, which implies 



C4- 



/i(C47) 
/o(C47)' 



(15) 



where /q and Ii are modified Bessel functions of the first 
kind. At the tetratic-nematic transition, we also have 
d^F/dC2 = 0, evaluated at C2 = 0, which implies 



7K 



7k/i(C47) 
/o(C47) 



(16) 



These two equations implicitly determine the second- 
order tetratic-nematic transition line shown in Fig. [31 
To find the tricritical point A, we expand the free en- 
ergy in powers of C2, for C4 satisfying Eq. psp . and we 
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require that the coefRcients of C| and C| both vanish. 
As a result, the tricritical point A occurs at 7 = 2.2496, 
K = 0.6116 and C4 = 0.4535. 

The first-order transition lines in the phase diagram 
cannot be calculated analytically; instead they are deter- 
mined by numerical minimization of the free energy. The 
triple point B occurs where the first-order transition lines 
intersect the second-order isotropic-tetratic transition of 
Eq. p^ . This point is found numerically at 7 = 2 and 
K = 0.79. 

Point C is the intersection of the extrapolated second- 
order transitions of Eqs. and p^ . which occurs at 
7 = 2 and k = 1. It represents the highest value of the 
symmetry breaking k where the tetratic phase can even 
be metastable, beyond the the triple point B where it 
ceases to be a stable phase. 



IV. MONTE CARLO SIMULATIONS 

So far, the calculations presented in this paper have all 
used mean-field theory. Of course, mean-field theory is an 



approximation, which tends to exaggerate the tendency 
toward ordered phases. In order to assess the validity of 
mean-field theory, we perform Monte Carlo simulations 
for a lattice model of the same system. In this lattice 
model, we use the Hamiltonian 

H = -J^{KCO8[2{0i - 6^)] + cos[A{0, ~ 0j)]) , (17) 

summed over nearest-neighbor sites i and j on a 2D tri- 
angular lattice, as shown in Fig. 01 This lattice Hamilto- 
nian corresponds to the model presented in the previous 
sections if we take the parameter 7 = 6J/{kBT), because 
each lattice site interacts with six nearest neighbors. 

We simulate this model on a lattice of size 100 x 100 
with periodic boundary conditions, using the standard 
Metropolis algorithm fllll. On each lattice site, the spin 
is described by an orientation angle 6. In each trial Monte 
Carlo step, a spin is chosen randomly, its orientation is 
changed slightly, and the resulting change in the energy 
AE is calculated. If energy decreases, the change is def- 
initely accepted. If not, the change is accepted with a 
probability of exp [—AE/{kBT)]. Usually, for a constant 
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temperature, each Monte Carlo cycle of the simulation 
consists of 10000 trial steps, and 50000 cycles are used for 
each temperature. However, near phase transitions, espe- 
cially near first-order transitions, additional Monte Carlo 
cycles are used to eliminate metastable states. The phase 
diagram is calculated by cooling the system from high 
temperature with decreasing the temperature in steps of 
0.01, or steps of 0.005 near phase transitions. During the 
last half of the simulation cycles, the order parameters 
are calculated and time-averaged. 

To calculate the nematic order parameter C2, we use 
the 2D nematic order tensor 

Qap = 2 {{uanp) - (naUp).^^) , (18) 

averaged over all lattice sites. Here, n ~ (cos 0, sin 0) is 
the unit vector representing each spin, and {nan/s)^^^ = 
^Sai3 is the average in the isotropic phase. The positive 
eigenvalue of this tensor is C2 . 

For the tetratic order parameters C4, we use the gener- 
alized tensor method of Zheng and Palffy-Muhoray [l^] ■ 
We consider the fourth-order tetratic order tensor 

= 4 [{naTipn^ns) - {naUpn^ns).^^^ , (19) 

averaged over all lattice sites. Here, we are subtract- 
ing off the isotropic average {nanpn^ns).^^^ = ^{Saf^S^s + 
Sa-fS/ss + SasSp^y). To Calculate the eigenvalues, we un- 
fold this fourth-order tensor into a second-order ten- 
sor (4x4 matrix), which we diagonalize using stan- 
dard methods. The four eigenvalues are 0, — C4, 
i (C4 - (16C| + C|)i/2), and i (C4 + (16C| + C|)i/2). 
(In the tetratic phase, with C2 — 0, they reduce to 0, 
— C4, +C4, and 0.) Thus, using the previously calculated 
value of C'2, we can extract C4. 

Figure [5] shows the simulation results for the order pa- 
rameters C2 and C4 as functions of 7, for several values of 
At. These results are quite simular to the numerical mean- 
field results of Fig. [21 although the quantitative values 
of 7, K, and the order parameters are somewhat differ- 
ent. For a small symmetry breaking k — 0.3, there are 
two second-order phase transitions. The order param- 
eter C4 increases continuously at the high-temperature 
isotropic-tetratic transition, and C'2 increases continu- 
ously at the lower-temperature tetratic-nematic transi- 
tion. For a slightly larger value of k = 0.5, the tetratic- 
nematic transition becomes first-order, with an appar- 
ently discontinuous increase in C2 (within the precision 
of the simulation). For k = 1, the intermediate tetratic 
phase disappears, and there is just a single first-order 
isotropic-nematic transition, with apparently discontinu- 
ous increases in both C2 and C4. Finally, for the largest 
value K — 3, the isotropic-nematic transition becomes 
second-order, with continuous increases in both C2 and 
C4. 

The simulation results are summarized in the phase 
diagram of Fig. [6] This phase diagram shows a high- 
temperature isotropic phase, an intermediate tetratic 
phase, and a low-temperature nematic phase. The tem- 
perature range of the tetratic phase is very large when 
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FIG. 6: Simulation results for the phase diagram in terms 
of 7 (inverse temperature) and k (two-fold distortion in the 
interaction). The triple point is at approximately 7 = 3.2 and 
K = 0.60. 

the symmetry breaking k is small, then it decreases as 
K increases, and eventually vanishes at the triple point, 
which is approximately given by 7 = 3.2 and k = 0.60. 
Compared with the mean-field phase diagram of Fig. [31 
the simulation results show the transitions at lower tem- 
perature (higher 7) than in mean-field theory. This dif- 
ference is reasonable because mean-field theory always 
exaggerates the tendency toward ordered phases. 

V. CONCLUSIONS 

In this paper, we propose a model for the statistical 
mechanics of particles with almost-four-fold symmetry. 
In contrast with earlier work on particles with a hard-core 
excluded- volume interaction, we consider particles with a 
soft interaction, analogous to Maier-Saupe theory of ne- 
matic liquid crystals. We investigate this model through 
two complementary techniques, mean-field calculations 
and Monte Carlo simulations. Both of these techniques 
predict a phase diagram with a low-temperature ne- 
matic phase, an intermediate tetratic phase, and a high- 
temperature isotropic phase. They make consistent pre- 
dictions for the order of the transitions and the tempera- 
ture dependence of the order parameters, although they 
do not agree in all quantitative details. 

The main result of this study is that the tetratic phase 
can exist up to a surprisingly high value of the symmetry 
breaking k in the microscopic interaction. We find the 
maximum k — 0.79 in mean-field theory, or 0.60 in Monte 
Carlo simulations. Even taking the lower Monte Carlo 
value, this implies that the interaction in the parallel 
direction (1 -I- k) can be about four times larger than the 
interaction in the perpendicular direction (1 — k). Hence, 
the tetratic phase can form even for particles with quite a 
substantial two-fold component in the interaction energy, 
i.e. for fairly rod-like particles. 
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It is interesting to note that our phase diagram is quite 
similar to the phase diagram found by density-functional 
theory for hard rectangles; see Fig. 3 of Ref. In that 
theory, the phase diagram shows isotropic, tetratic, and 
nematic phases, and the tetratic phase can exist for rect- 
angles with aspect ratio of up to 2.21:1. That theory 
shows the same arrangement of the phases, and even 
the same first- and second-order phase transitions, with 
tricritical points on the isotropic-nematic and tetratic- 
nematic transition lines. This phase diagram seems to 
be a generic feature of particles with four-fold symmetry 
broken down to two-fold. Thus, we can expect to see 



tetratic phases in 2D experiments and simulations, even 
if the particles are moderately extended. 
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